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Abstract 

We examine the quantum correlations of spin pairs in the ground state of finite XY chains in 
a transverse field, by evaluating the quantum discord as well as other related entropic measures 
of quantum correlations. A brief review of the latter, based on generalized entropic forms, is also 
included. It is shown that parity effects are of crucial importance for describing the behavior of 
these measures below the critical field. It is also shown that these measures reach full range in 
the immediate vicinity of the factorizing field, where they become independent of separation and 
coupling range. Analytical and numerical results for the quantum discord, the geometric discord 
and other measures in spin chains with nearest neighbor coupling and in fully connected spin arrays 
are also provided. 
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I. INTRODUCTION 



The last decades have witnessed the great progress experienced by the interdisciplinary 
field of quantum information science jl-3], which began with the recognition of the potential 
of quantum systems and quantum correlations for information processing tasks. While it 
is well known that quantum entanglement is essential for quantum teleportation |4j, super- 



dense coding, [5[ and also for achieving exponential speed-up in pure state based quantum 
algorithms [6|, the mixed state based quantum algorithm of Knill and Laflamme p| showed 
that such speedup could in principle be achieved in this case without a substantial presence 
of entanglement 8j . This has oriented the attention to alternative definitions and measures 
of quantum correlations for mixed states, like the quantum discord jsl, Q . While coinciding 
with the entanglement entropy for pure states, the quantum discord differs essentially from 
the entanglement of formation in the case of mixed states, being non-zero in most separable 
mixed states and vanishing just for states which are strictly classically correlated at least 
with respect to one of the constituents, i.e., diagonal in a standard or conditional product 



basis [9]. The result of Ref. ll| showing the existence of a finite discord between the eon1 vol 
qubit and the remaining qubits in the circuit of Ref. jl], unleashed a great interest on tins 



measure and several investigations on its fundamental properties |12Ml5|, on its evaluation 
on spin chains and specific states 16M23| as well as on related measures 24j-l29j]. have been 
recently made (see Ref. 30] for a recent review). Distinct quantum capabilities of states 



with non-zero discord have also been recently investigated 

Our aim here is to describe the remarkable behavior of the quantum discord and of other 
related entropic measures of quantum correlations, in the exact ground state of finite XY 
chains in a transverse field 20j. We first provide in sections 2-5 a brief review of the quantum 



discord and of the generalized entropic measures of quantum correlations discussed in Refs. 
28 . 2^] . The latter comprise as particular cases the one-way information deficit 31, 35 1 
and the geometric measure of discord of ref. [27], embedding them in a unified formalism 



based on majorization 



36 



371 ] and general entropic forms 



38 



While their basic features 



are similar to those of the quantum discord, the possibility of using simple entropic forms 
permits an easier evaluation, allowing for analytical expressions in some cases, as occurs 
with the geometric discord of general two qubit states 271 ] . 

We then use these measures to investigate, in sections 6-8, the quantum correlations of 



spin pairs in the exa ct g round state of finite XY chains in a transverse field. We review the 
main results of Ref. on the behavior of the quantum discord in these chains and also add 
new results concerning the behavior of the geometric discord and other related measures in 
such chains. The exact ground state of a finite XY chain in a transverse field has a definite 
spin parity and this fact will be seen to deeply affect the discord and the previous measures 
for fields lower than the critical field B c . We will show that the essential results in this sector 
can be interpreted in terms of the discord of mixtures of aligned pairs. 



Moreover, these chains can exhibit a factorizing field B s [39H47I]. where they have a 
completely separable ground state. For transverse fields, such eigenstate actually breaks 
the previous parity symmetry and is hence degenerate, coinciding B s in a finite chain with 
the last crossing of the two lowest opposite parity levels 43j. A related remarkable effect 
is that in the immediate vicinity of B s , pairwise entanglement, though weak, reaches full 



range 42|, |43], regardless of the coupling range. Here we will show that the quantum discord 



as well as the entropic measures of quantum correlations also reach full range at this point, 
exhibiting universal features such as being independent of separation and coupling range 



201 ] . Moreover, the value reached by them at this point is non-negligible and does not 



decrease with size, in contrast with the pairwise entanglement, since these measures are not 
restricted by the monogamy property [48j which affects the latter (limiting the concurrence 
(49! to order n _1 in an n spin chain if all pairs are equally entangled). Consequently, the 
behavior of these measures with the applied field and separation will deviate significantly 
from that of the concurrence or entanglement of formation for \B\ < B s . Conclusions are 
finally discussed in section 9. 



II. QUANTUM DISCORD 

The quantum discord was originally defined {9], [uj] as the difference between two distinct 
quantum versions of the mutual information, or equivalently, the conditional entropy. For 
a classical bipartite system A + B described by a joint probability distribution pij = p(A = 
i,B = j), the conditional entropy is defined as the average lack of information about A 
when the value of B is known: S(A\B) = j pf ''S \A\B = j), where pf = J2iPij is the 
probability of outcome j in B and S(A\B = j) = — J2iPi/j ^°SPi/j is the Shannon entropy 
of the conditional distribution pi/j = Pij/pf ■ It is a non-negative quantity, and can also be 
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expressed in terms of the joint entropy S(A, B) = — J2i j Pij ^°&Pij an d the marginal entropy 
S(B) = -J2jPj l °SPf as S(A\B) = S(A,B) - S(B). Positivity of S(A\B) then implies 
S(A,B) > S(B) (and hence S(A,B) > S(A)) for any classical system. 

The last expression for S(A\B) allows a direct quantum generalization, namely 

S(A\B) = S(p AB ) - S(p B ) , (1) 

where S(p) = — Trplog 2 p is now the von Neumann entropy and pab the system density 
matrix, with p B = Tr a Pab the reduced state of subsystem B. It is well known, however, 



that Eq. ([I]) can be negative 36|, being for instance negative in any entangled pure state: 
If p \ B = PAB) S(pab) = and S(A\B) = —E(A, B), where E(A, B) = S(A) = S(B) is the 
entanglement entropy [50]. The positivity of Eq. ([rj provides in fact a basic separability 



criterion for general mixed states 5lJ: pab separable =>- S(A\B) > 0, a criterion which can 



actually be extended to more general entropies 52 



541 ] . We recall that pab is separable if it 



can be written as a convex combination of product states, i.e. pab = ^2 a QaP A ® P%i "with 



Qa > 0, Y.aQa = 1 [55 1- 



A second quantum version of the conditional entropy, closer in spirit to the first classical 
expression, can be defined [s| on the basis of a complete local projective measurement M B 
on system B (von Neumann measurement), determined by one dimensional orthogonal local 
projectors Hj = \]b){3b\- The conditional entropy after such measurement is 

S Mb (A\B) = Y.Pf S (pA/i) = S(p'ab) - S(p' B ) , (2) 

3 

where pf = Trp^Il^, with IT^ = I A <8> II.,-, is the probability of outcome j, pA/j = 
Tr B (pA B Hf)/pf is the reduced state of A after such outcome and 

p'ab = E P?PW ® U i = J2 U fPAB^f , (3) 

3 j 

is the average joint state after such measurement, with p' B = Tr a p'a B = YljPf^-j- Eq. 
(J2J) represents the average lack of information about A after the measurement M B in B 
is performed and is clearly non- negative, in contrast with Eq. (CQ). For a classical system, 
both quantities (DO)-® are, however, equivalent. We also mention that for general local 
measurements, defined by a set of positive operators Ej Ej = I), the first expression in 
(f2]) is to be used for S Mb {A\B) (with EL, Ej). 



The quantum discord 9Mlll| can then be defined as the minimum difference between Eqs. 
© and flU): 

D B (p AB ) = Mm[S MB (A\B)] - S(A\B) , (4) 

Mb 

where the minimization is over all local measurements M B - Due to the concavity of the 
von Neumann conditional entropy (pQ) with respect to pab, [36( Eq. @ is non- negative 
9|, vanishing just for classically correlated states with respect to B, i.e. states which are 
already of the general form ([3]) (a particular case of separable state) and which remain 
then unchanged under a specific unread local measurement. Such states are diagonal in a 
"conditional" product basis = \ij A ) <8> Ub)}, with \ij A ) the eigenstates of PA/j- 

Eq. (jl]) is then non-zero not only in entangled states but also in separable states not of 
the form ([3]), i.e. those which involve convex mixtures of non-commuting product states, for 



which the entanglement of formation [56( vanishes. It is then a measure of all quantum-like 
correlations between A and B. The distinction with entanglement arises nonetheless just 
for mixed states: For pure states p AB = pab, S(pab) = and S(p' AB ) = S(p' B ) for any von 
Neumann measurement, reducing the quantum discord exactly to the entanglement entropy: 
D A = D B = E(A, B). 

Eq. (J3J) can of course be also understood [h]] as the minimum difference between the 
quantum mutual information [36( I(A : B) = S(A) — S(A\B) = S(pa) + S(p B ) — S(pab), 
which measures all correlations between A and B (I (A : B) > 0, with I (A : B) = 
if and only if pab — Pa 8> Pb) and the "classical" mutual information Im b {A : B) = 
S(A) — Sm b {A\B), which measures the correlations after the local measurement Mb- 



III. GENERALIZED ENTROPIC MEASURES OF QUANTUM CORRELATIONS 



Let us now discuss an alternative approach for measuring quantum correlations [281 ]. 
which allows the direct use of more general entropic forms. We consider a complete local 
projective measurement Mb (Von Neumann type measurement) on part B of a bipartite 
system initially in a state pab, such that the post-measurement state is given by Eq. (j3J) if 
the result is unread. A fundamental property satisfied by the state (j3J) is the majorization 
relation [l|, 0, 37] 

Pab ~< Pab , (5) 



where p' -< p means, for normalized mixed states p, p' of the same dimension n, 



3=1 3=1 



l,...,n- 1 



Here p^-, pj denote, respectively, the eigenvalues of p' and p sorted in decreasing order (p^ > 0, 
J2jPj = !)■ Eq. (EJ) implies that p' AB is always more mixed than p A #: If Eq. ([HJ) holds, p AS 



can be written as a convex combination o: 
q a >0,J2 a q a = l and W a U a = I 13 , 



37|. 



unitaries of p AB : P^b = E Q ^^Pab^, with 



Eq. (JSJ) not only implies that S(p' AB ) > S(pab) for the von Neumann entropy, but also 



Sf(p' AB ) > Sf(pAB) 



(6) 



for any entropy of the form 38 ] 



S f (p)=Trf(p), (7) 

where / : [0, 1] — > 3? is a smooth strictly concave function satisfying /(0) = /(l) = 0. 
As in the von Neumann case, recovered for f(p) = — plogp, these entropies also satisfy 
Sf(p) > 0, with Sf(p) = if and only if p is a pure state (p 2 = p), and S/(p) maximum for 
the maximally mixed state p = I n /n. Hence, Eq. (j3j) implies a strict disorder increase by 
measurement which cannot be fully captured by considering just a single choice of entropy 
(S(p') > S(p) does not imply p' -< p). More generally, Eq. (jSJ) actually implies F(p') > F(p) 
for any Schur concave function F of p[37(]. Nonetheless, entropies of the form (13) are sufficient 
to characterize Eq. (jHJ), in the sense that if Eq. (jHJ) holds for a// such Sf, then p' -< o p4 \. 
We may now consider the generalized information loss due to such measurement 2s| - 



I?* = S/Wju*) ~ Sf(pAB) 



(8) 



which is always non-negative due to Eqs. (jS])-©? vanishing only if p' AB = pab due to the 
strict concavity of /. Eq. (jSj) is a measure of the information contained in the off-diagonal 
elements (ij\pAB\i'j') {j ^ f) of the original state, lost in the measurement. The minimum 



of among all complete local measurements 



28|, 



If ( P ab) = Mm S f (p' AB )-S f (pab) 

Mb 



(9) 



provides then a measure of the quantum correlations between A and B present in the original 
state and destroyed by the local measurement in B: If > 0, vanishing, as the quantum 



discord (j3J), only if p AB is already of the form ([3]), i.e., only if it is diagonal in a standard or 
conditional product basis. 

Again, in the case of a pure state p AB = \^ ab){^ ab\, it can be shown |28| that Eq. 
((9]) reduces to the generalized entanglement entropy: If = if — Ef(A,B) if p\ B = pab, 
where Ef(A,B) = Sf(p A ) = Sf(ps)- The minimizing measurement in this case is the local 



Schmidt basis for \V AB ): M B = {\k B }(k B \} if \V AB ) = Ek y/Pk\ k > ® \ k ) 



28 



In the case of the von Neumann entropy (Sf(p) = S(p)), Eq. <Q becomes the one-way 



information deficit |3lL l35|. which coincides with the different version of discord given in the 



last entry of Ref. [9[ (and denoted as thermal discord in [30|]). It can be rewritten in this 



57j S(p\\p') = — Trp(logp' — logp) (a non-negative 



case in terms of the relative entropy [36 . 
quantity) as 

I b ( P ab) = Mm S(p' AB ) - S(pab) = Mm S(p AB \ \p'ab) , 



(10) 



M B " M B 

where we have used the fact that the diagonal elements of pab and p' AB in the basis where 
the latter is diagonal are obviously coincident. We also note that for these measurements, 
the quantum discord (@| can be expressed as D B = Min M B [l M B (pab) —I Mb {pb)]i coinciding 
then with I B when the minimizing measurements in ffTUl) and (jl]) are the same and such 
that p' B = Pb- 

In the case of the linear entropy S^p) = 1 — Trp 2 , obtained for /(p) = p(l — p) (i.e., 
using the linear approximation logp — >• p — I in — plogp), Eq. flS]) becomes 



h (Pab) = Mm Tr(p AB - p AB ) = Mm \\p' AB ~ Pab\ 

Mb Mb 



(11) 



where ||0|| 2 = Tr 0^0 is the squared Hilbert Schmidt norm. This quantity becomes then 
equivalent to the geometric measure of discord introduced in Ref. j^ . The latter is defined as 
the last expression in Eq. ( II ip with minimization over all states diagonal in a product basis, 
but the minimum corresponds to a state of the form ([3]) [28]. For pure states, I B becomes 
proportional to the squared concurrence C\ p , 49] as for pure states C\ B is proportional to 



the linear entropy of any of the subsystems 58 



Finally, in the case of the Tsallis entropy [59J S q (p) = (1 — Trp 9 )/(g — 1), q > 0, which 
corresponds to f(p) = (p — p q )/{q — 1), Eq. (j9]) becomes [29] 



IqipAB) = Mm S q (p' AB ) - S q {p AB ) oc MinTr (p q AB - p' q AB 



M 



Mr 



(12) 
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with I q reducing to the one way information deficit for q — > 1 (as S q (p) —> S(p)), and to 
the geometric discord for q = 2. This entropy allows then a simple continuous shift between 
different measures. 

When considering qubit systems, we will normalize entropies such that Sf(p) = 1 for 
a maximally mixed two-qubit state p (i.e., 2/(1/2) = 1), implying that all if will take 
the value 1 in a maximally entangled two-qubit state (Bell state). This implies setting 
log = log 2 in the von Neumann entropy, S^p) = 2(1 — Tr p 2 ) in the linear case (such that 
If = 2Min A f B ||p^ B - Pab\\ 2 ) and S q (p) = (1 - Trp 9 )/(1 - 2 1 - q ) in the Tsallis case. 



IV. GENERAL STATIONARY CONDITIONS FOR THE LEAST DISTURBING 
MEASUREMENT 

The stationary condition 5I^ Ib (pab) = for the quantity (|S]), obtained by considering a 
general variation S\js) = (e lSflB — I)\jb) ~ iShljs) of the local measurement basis, where 



Kb is an hermitian local operator, reads 



29|] 



Tr A [f'(p' AB ),p AB } = 0, (13) 

i- e -> ^L,i[f'(Pj)(ijj\PAB\ijk) - f'(pl)(ikj\PAB\ikk)] = 0, where /' denotes the derivative of / 
and {ijj\pAB\i'jj) = $a'P % j- In the case of the geometric discord, f'(p) oc / — 2p and Eq. (|T3|) 
reduces to Tr^ [p'abiPab] = 0. In the case of the quantum discord (j4j), Eq. ([TBI should be 



replaced for these measurements by 29] 



Tr A [f (Pab),Pab] - [f(p' B ),PB} = 0, (14) 

with f(p) = — plogp, due to the extra local term. 

Eqs. (fTBl - fTHl) allow us to identify the stationary measurements, from which the one 
providing the absolute minimum of I^ Ib (least disturbing measurement) is to be selected. For 
instance, if there is a standard product basis where (ij\pAB\if) = fyj'Pj an d (ij\pAB\i'j) = 
dwPj, such that the only off-diagonal elements are {ij\pAB\i'j') with i ^ i' and j ^ j', 
a measurement in the basis {\j B )} is clearly stationary for all If, as Eq. (jTBl is trivially 



291 ] . It will also be stationary for the 



satisfied, leading to a universal stationary point 
quantum discord. 

An example of such basis is the Schmidt basis for a pure state, \^ab) = 2^k=i \/Pk\kk), 
with \kk) = \k A ) <S> \ks) and n s the Schmidt rank, since (kl\pAB\k'V) = ^ki^k<v s/PkPk 1 for 
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Pab = \^ab) (&ab\- The same holds for a mixture of a pure state with the maximally mixed 

state, 

p AB (x) =x\^ AB )(q AB \ + —^I n , (15) 
where n = n A n B and x G [0,1]. The Schmidt basis provides in fact the actual minimum 



of lf(x) = If (pab(x)) V x G [0,1], as shown in Ref. [5 



. This implies the existence in 



this case of a universal least disturbing measurement, and of a concomitant least mixed post 
measurement state, such that p' AB majorizes any other post-measurement state. We can 
then obtain a closed evaluation of If for this case V Sf, [28] which shows some of its main 
features: 

B = £ *(^-l) + l _ »(n-l) + l _ _ 1-* 

k=l 

If n s > 1, it can be shown that lf{x) > for x > 0, being a strictly increasing function 
of x for x G [0, 1] if /(p) is strictly concave. Moreover, a series expansion for small x leads 
to If(x) ~ ax 2 (l — X/fc=iPfc)' wnere a = ~il"(^/ n )/^ > 0, indicating a universal quadratic 
increase with increasing x if /"(1/n) 7^ [28]. This behavior is then similar to that of the 
quantum discord [9| and quite distinct from that of the entanglement of formation, which 
requires a finite threshold value of x for acquiring a non-zero value. 



V. THE TWO QUBIT CASE 

Let us now examine the particular case of a two qubit system. A general two-qubit state 
can be written as 

pAB = ^{I + r A -o- A + r B -er B + cr A J(T B ), (17) 

where I = I2 <E> I2 denotes the identity, a A = a <g) I2 and a B = I2 <S> c. Due to the 
orthogonality of the Pauli matrices, we have r A = (cr A ), r B = (cr B ) and J = (cr A er B ), i.e., 
J m ' = (ctah<7bh')i where p, p! = x,y,z and (O) = Tr p AB O. 

A general local projective measurement in this system is just a spin measurement along 
a unit vector k, and is represented by the orthogonal projectors |(/± fc • cr). Therefore, the 
most general post-measurement state (j3J) reads 

Pab = \{I + r A - <t a + (r B ■ k)k ■ a B + {(T A Jk){k ■ a B )) , (18) 

9 



and corresponds in matrix notation (setting r and k as column vectors) to rg kkfrs and 
J — > Jkk 1 . The general stationary condition f lT3|) can be shown to lead to the equation j^j 



air B + a 2 J t r A + a 3 J l Jk = Xk , (19) 

i.e., k x (airs + cx2J t fA + cis^Jk) = 0, which determines the possible values of the mini- 
mizing measurement direction k. Here (ai,a 2 ,a 3 ) = \ ^2 /'(p")^ 1 ' 1 ''/^)^/^)) with 

i/,i/=±i 

=j(l + z/r B -fc + z/'Aj,) the eigenvalues of ffTSl) and = + ^Jfe|. In the case of the 



quantum discord (j3]), the additional local term leads to the modified equation 



(a 1 -7])r B + a 2 J t r A + a 3 J t Jk = Xk , (20) 

where here f(p) = -plogp and rj = |^=±i^ 7' '{Pv) = \^>g{p-/p + ), wi th Pv = J2„>Pv = 
|(1 + vtb ■ k) the eigenvalues of p' B . A different approach was provided in Ref. [2^. 

General analytic solutions of these equations can be obtained in a few cases. For instance, 
a closed evaluation of If for any Sf is directly feasible for any two-qubit state with maximally 
mixed marginals, i.e. 

Pab = \( J + VaJvb) , (21) 
for which Eq. (Tl9|) reduces to J l Jk = Xk V If, indicating that k should be an eigenvector of 



J 1 J. Moreover, it can be shown 29[ that the minimum corresponds to k directed along the 
eigenvector with the largest eigenvalue of J* J V Sf (universal least disturbing measurement), 
such that the post-measurement state (l3l)-(|T8l) conserves the largest component: By suitable 
local rotations, <taJ(t b can be written as Y^a=xyz J^a^b^ where are the eigenvalues 
of J* J (the same as those of J J 1 ), and the least disturbing measurement leads then to 
p'ab = \{I + JfL (J Ajj. (J Bii, where Jp = Max[J s , J y , J z ]. The final result for If (obviously 
identical to 1^ for this state) is 29] 



If {Pab) = 2f(p^) + 2f(^±2i) - f( Pl ) - /(pa) - f(p 3 ) - /(p 4 ) , (22) 



l + J Z ±(J X -Jy) 
4 ' 



where ( P i, P 2, P 3, Pa) are the eigenvalues of pab sorted in decreasing order (p 12 
p 3A = i-J'^+Jy) if \ j z \ > \j x \ > \ j y \ an d j z > 0, J x > 0). It is verified that If = only 
if p\ = P2 and p 3 = p 4 , in which case pab = Pab 1S a classically correlated state (29] . 
In the von Neumman case, Eq. (122]) is just the quantum discord for this state, as in this 
case it coincides with the one-way information deficit (p' B = p B are maximally mixed). In 

10 



the case of the linear entropy, Eq. (122]) yields the geometric discord and reduces to [29] 

IHpab) = (Pi - P2) 2 + (p 3 - P4) 2 . 

The linear entropy case fill I) is obviously the most simple to evaluate, and in this sense 
the most convenient. A full analytic evaluation for a general two-qubit state was achieved 
in Ref. 



271 ]. Since Tr cr^ay = 28^, one easily obtains in this case 



I2 = S2(p' AB ) - S 2 (p AB ) = ~(trM 2 - k t M 2 k) , M 2 



fBr B 



R + J f J. 



(23) 



where ||J|| 2 = trJ*J, |r| 2 — r ■ r and M 2 is a positive semidefinite symmetric matrix. The 
minimum of J 2 is then obtained when k is directed along the eigenvector associated with 
the maximum eigenvalue Ai of M 2 , leading to 27] 

1 



if = Mm/ 2 fc 



-(trM 2 -Ai) 



(24) 



It is easily seen that Eq. f fl9|) reduces in this case to the eigenvalue equation M 2 k = Xk, 
such that the stationary directions are those of the eigenvectors of M 2 . 

Similarly, the q = 3 case in the Tsallis entropy, S 3 (p) ex (1 — Trp 3 ), can also be fully 
worked out analytically 29]. We obtain 



rk 

M 3 



S 3 (pab) = ^(tr M 3 - 2det J - k l M z k) 



SMb 

r B r B + J* J + r B r A J + J l r A r B 



(25) 
(26) 



where M 3 is again a positive semidefinite symmetric matrix, with trM 3 = |t*_b| + ||^|| + 
2r A Jr B - Its minimum corresponds then to k along the eigenvector with the maximum 



eigenvalue Ai of M 3 , which leads to 



29] 



J 3 B = Min / 3 fe = -(tr M 3 - 2det J - Ai) . 



(27) 



It is again verified that Eq. ([19]) leads here to the same eigenvalue equation M 3 fc = Xk, 
as (ai,a 2 ,a 3 ) = {r B + r A J,r B k,l). In the case of the state (12T1) . we obtain [29| /g 8 = 
(pi - p 2 ) 2 (pi + p 2 ) + (p 3 - p 4 ) 2 (P3 + P4)- 

It should be stressed that for the case of two qubits, these two entropies, S 2 and S 3 , 
lead to the same entanglement monotone [6Q], since for an arbitrary single qubit state they 



become identical: 



28 



29| S 2 ( PA ) = S 3 (p A ) = 1 - M 2 for p A = \{I 2 + r A 



Both 



quantities I2 and 1% reduce then to the standard squared concurrence |49| C\ B in the case 
of a pure two-qubit state. 
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VI. THE CASE OF A MIXTURE OF TWO ALIGNED STATES 



We are now in a position to examine the important case of a mixture of two aligned spin 



1/2 states 20j, which will allow us to understand the behavior of the quantum discord of 



spin pairs in finite ferromagnetic XY chains, particularly in the vicinity of the transverse 



factorizing field 



20 



39 



47J. We consider the bipartite state 



Pab(9) = \{\ee)(M\ + \-0- 9)(-9 - 9\) (28) 
= -(I + cos 9(o Az + a Bz ) + cos 2 9 o Az o Bz + sin 2 9 a Ax a Bx ) , (29) 

where \9) = e~ ie °~ y / 2 \0) denotes the state with its spin aligned along k = (sm9, 0, cos (9), and 
( 1291) corresponds to the spin 1/2 case, where \9) = cos ||0) + sin ||1) and we have used the 



22j, i.e., states that commute with 



notation of Eq. ()T7|) . It is a particular case of X state 
the S z parity P z = - e i7r ( CT ^ +CT ^)/ 2 [29j. 

The state (1281) arises, for instance, as the reduced state of any pair in the n-qubit pure 

states 

. n , \99...9)±\-9...-9) 
y/2(l ± (-9\0)») 

if the complementary overlap (— 9\9) n ~ 2 ((—9\9) = cos 9 for spin 1/2) can be neglected (i.e., 
n large and 9 not too small). As will be seen in the next sections, the states ( 1301) are the 
actual exact ground states of such chains in the immediate vicinity of the factorizing field. 



The state (128!) is clearly separable, i.e., a convex combination of product states 55], 
but is classically correlated, i.e., diagonal in a product basis, just for 9 = or 9 = tt/2. 
Accordingly, both the quantum discord and all measures including the geometric discord, 
will be non-zero just for 9 G (0,7r/2). As seen in Fig. [TJ they all exhibit similar qualitative 
features, although significant differences concerning the minimizing measurement arise. Due 
the symmetry of the state, it is apparent that D A = D B = D and I A = If = If V 9. 

It is rapidly seen from Eq. (1191) that for this state (as well as any other X state), spin 
measurements along x, y or z are stationary, for both the quantum discord and all measures 
If {2^]. In the case of the quantum discord, the minimizing measurement for this state 
which is of rank 2, and hence minimized through a standard von Neumann measurement 
3(3]) is in fact along x V 9 E (0,7r/2), in which case the eigenvalues of p' AB become p^' = 
|(1 + vW cos 2 9 + sin 4 9), being twofold degenerate. The final result for the quantum discord 
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71/4 7C/2 



FIG. 1. Quantum correlation measures in the mixture of aligned states ()28|) : The quantum discord 
D, the geometric discord J 2 and the "cubic" discord I3, as a function of the angle 0. Normalization 
is such that all measures take the value 1 in a maximally entangled two-qubit state. Due to the 
symmetry of the state, D = D A = D B and If = 1^ = if V /. 



can then be expressed as |20j 



1 + vJ 1 - \ sin 2 20 1 +//cos 2 6> l + z/cos0 



u=±l 



where f(p) = —plog 2 p- It is maximum at 9 ~ 1.157r/4. For 9 ~ 0, D vanishes quadratically 
(D oc 9 2 ) whereas for 9 ->■ tt/2, D cx (f - ^) 2 (- log 2 (| - ^) 2 + c). 

On the other hand, the geometric discord (1111) and the "cubic" discord (q = 3 in Eq. 



dHJ) can be directly evaluated using Eqs. (l2i|) -([27Il . We obtain [29| 



i sin 4 < 0, 



2 



c2 



I cos 2 + cos 4 > c2 

I sin 4 < 0, 

i(cos 2 + 3cos 4 0) 0>0, 



c3 



(32) 



(33) 



c3 



where cos 2 c2 = 1/3 (0 c2 w 0.61tt/2) and cos 2 c3 = (v 7 !^ - 3)/4 (0 c3 w 0.64vr/2), the 
minimizing measurement direction changing abruptly from z to x at 9 = 9 C as increases, in 
contrast with the quantum discord. Both J 2 and ^3 exhibit therefore a cusp like maximum at 
= C , as seen in Fig. [H It is also seen from Eqs. ( I32l -(l33"l) that these quantities vanish as 4 
for — y 0, whereas for — > 7r/2 they vanish quadratically (oc (n/2 — 9) 2 ). For completeness, 
it should be also mentioned that the behavior of the least disturbing local measurement for 
this state depends actually on the choice of entropy. For instance, in the von Neumann case 
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(where Eq. (jUJ) becomes the one-way information deficit ([TO]) ), we obtain instead a smoothed 
z — > x transition for the minimizing measurement direction, which evolves continuously from 
z to x in a small intermediate interval 29]. 



VII. XYZ SPIN CHAINS AND TRANSVERSE FACTORIZING FIELD 



Let us now use the previous measures and results to analyze the quantum correlations 
between spin pairs in a chain of spins Sj. We will consider finite chains with XYZ couplings 
of arbitrary range immersed in a transverse magnetic field, not necessarily uniform, such 
that the Hamiltonian reads 



H = J2B % s lz -\ J2 E J ? s ^>> (34) 

i /i=x,y,z i,j 

where Sj M denote the components of the local spin (assumed dimensionless) . We first 
remark that the Hamiltonian ( 1341) always commutes with the total S z parity or phase flip 



>i=l 



exp[— nr[Si 



Si)} 



(35) 



irrespective of the coupling range, anisotropy, geometry or dimension of the array. Hence, 
non-degenerate eigenstates will have a definite parity In fact, the ground state of finite 
chains will typically exhibit a series of parity transitions as the field increases from 0, before 
ending in an almost aligned state for sufficiently large fields. 

A related remarkable effect in these chains is the possibility of exhibiting a completely 
separable exact eigenstate at a factorizing field. The existence of a factorizing field was 
first discussed in Ref. j^ . and its properties together with the general conditions for its 



existence were recently analyzed in great detail by several authors 20|, I4CH47I] . At the 



transverse factorizing field, finite XYZ chains actual 
and degenerate parity breaking exact eigenstates [43 



y exhibit a pair of completely separable 



45[ , which can be ground states under 



quite general conditions. In such a case the transverse factorizing, field corresponds to a 



ground state parity transition (typically the last parity transition 43 



where the lowest 



energy levels of each parity subspace cross and enable the formation of such eigenstates. Let 
us notice that while these lowest levels become practically degenerate in a large chain for 
fiel ds \B\ < B c , they are not exactly degenerate in a finite chain, except at crossing points 



43 



45 



61 



62 1 
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Let us then first describe the general conditions for which a separable parity breaking 
state of the form 



|9> = |0i ...d n ) = ^xexph^s^lO,) 



(36) 



where Sj Z \0j) = —Sj\0j), can be an exact eigenstate of 
H\Q) = E\Q), it can be shown that such conditions are 



'. 3y inserting in the equation 



45| 



Jy — Jx COS ^ COS % + Jz Sil1 ^ Sm ®3 ' 

5* sin 6*4 = — ~^ij){Jx cos ^* sm % — ^ s i n ^ cos ^j) > 



(37) 
(38) 



which are valid for arbitrary spins Sj. They determine, for instance, the values of J l J and B l 
in terms of J l J , J l J , and 6*j. A careful engineering of couplings and fields can then always 
produce a chain with such eigenstate, for any chosen values of Q{. It is also apparent that 
this eigenstate is degenerate, since P Z \Q) = \ — 0) = ®"=i ex P[^j s jy\\^j) wm have the same 
energy (and differ from |0) if sin0j ^ for some j), indicating that these fields necessarily 
correspond to the crossing of two opposite parity levels. Each local state in the product (|36l) 



is a local coherent state. We also note that any state ®j =1 e _:(0J " Sj \0j) can be written, except 
for a normalization factor, in the form ( 1361) by allowing a complex 9j 45] . Eqs. (I3~7l)-(l38l) 
are then generally valid for such type of states. 

The second equation ( 138]) cancels the matrix elements of H between |0) and one spin 
excitations, and is then a "mean field- like" equation, i.e., that which arises when minimizing 
the average energy (0|if|0) with respect to the 0j, for fixed fields and couplings. The 
first equation ( )37|) ensures that the minimizing separable state is an exact eigenstate, by 
cancelling the residual matrix elements of H connecting |0) with the remaining states (two- 
spin excitations). It can be also shown 45j that in the ferromagnetic- type case 



\4 j \<Ji j v/../ 



(39) 



where all off-diagonal elements of H in the standard basis of s z eigenstates are real and 
negative, the state ( 1361) is necessarily a ground state if 9j G (0, 7r/2) V j, as the exact ground 
state must have (or can be chosen to have if degenerate) expansion coefficients of the same 
sign in this basis (different signs will not decrease (H)), and hence cannot be orthogonal to 



I©) mm 
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In particular, a uniform solution 9j = V j, leading to |G) = \6 ...6), is feasible if the 
coupling anisotropy 

\ TT-rJ- ( 4 °) 

is constant V i, j and non-negative 



45j. Of course, if x > 1 we can change it to x G (0, 1) 



by swapping x •<-)■ y through a rotation of tc round the z axis. In such a case, Eqs. ( |37|) -( l38|) 
lead to 



cos 2 = x, 



(41) 
(42) 



where Eq. ( 1421) holds for sin 8 0. 

Eqs. fj4T]) - fr42|) allow, for instance, the existence of a factorizing field for uniform first 
neighbor couplings J l J = J^Sj^i in a finite linear spin s chain if x 



Jy — Jz 



cyclic case (J* n = J^), where the factorizing field is completely uniform, 

B l = B s = 2s^x~(J x - J z ) , 



> 0, both in the 



(43) 



as well as in the open case (J^ n = 0), where B l = B s at inner sites but B 1 = B n = B s /2 at 
. A fully and equally connected spin s array with = 2J A1 /(n — 1) V i ^ j 



the borders 



(Lipkin model 63|, |64|) will also exhibit a uniform transverse factorizing field at B = B s if 



X > |43j, |65(. The ensuing state |0) can be ensured to be a ground state in all these cases 



if | Jy | < Jx (when \ G [0, 1]). Other possibilities, like solutions with alternating angles 45], 
can also be considered. 



VIII. QUANTUM CORRELATIONS IN THE DEFINITE PARITY GROUND 
STATES 

Let us now focus on finite spin chains which exhibit a separable parity breaking exact 
eigenstate |0) at the factorizing field B s . It will be of course degenerate with | — O) = P Z \Q). 
The important point is that the definite parity states 

|8 ±>= ' 9)±| - 9> . (44, 

V2(i ± <-e|e» 

i.e. the states ( |30|) in the uniform case, will also be exact ground states at B s . Moreover, 
since the exact ground state of a finite chain will actually be non-degenerate away from the 
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factorizing field (and the other crossing points), it will have a definite parity. Hence, the 
actual ground state side-limits at B = B s will be given by the definite parity states (jUJ) 
(rather than | ± 6)). A ground state parity transition > + will then take place as the 



field increases across 



43 



45] 



While the ground states of each parity sector become degenerate in the large n limit for 
fields \B\ < B s , in finite chains the degeneracy is lifted and the actual ground state will 
exhibit important correlations arising just from the definite parity effect. For instance, in 
the immediate vicinity of B s , the pairwise correlations, rather than vanish, will approach 
the values determined by the states (jSJ). They will then depend on 0j and 9j for a pair 
i,j, irrespective of the separation i — j. In the uniform case, such correlations will then be 
independent of the separation, since the states (jUJ) will be completely symmetric and will 
lead to a separation independent pair reduced state pij = Trj--|6±)(6±| (ij denotes the rest 
of the chain). Such state is given exactly by 20| 

\eo)(oe\ + \-9-9){-e-e\+e{\ee){-e-e\ + \-e-e)(ee\) 



(45) 



2(l + £ (00|-0-0)) 

where e — ± cos n_2 for s = 1/2. This parameter is then small for not too small n and 
In the e — > limit one recovers the mixture fl28l). 



The concurrence of the state (141)1) (a measure of its entanglement 49[) depends essentially 



on e and is therefore small, vanishing 



separable). It is given explicitly by 
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'or e — > (as in this limit the state (T4"5l) becomes 



43 



45| 



C 



\e\ sin 2 



(46) 



1 + e cos 2 ' 

which is parallel (antiparallel) |42] for e > (< 0). Its maximum value is 2/n (in agreement 
with the monogamy property 481]), reached for 9 — > in the negative parity case, where the 
state ( |30l) approaches an W-st&te |43l |. 

In contrast, we have seen that the quantum and geometric discord, as well as the other 
measures of quantum correlations (Q, do acquire finite and non- negligible values in the 
mixture (1281) (Fig. [1]), i.e., in the state (T4"5l) even for e — Y 0, entailing simultaneous and 
coincident finite values for all pairs in the state ( 130]) 201 ] . This implies in turn infinite range 
of pairwise quantum correlations, as measured by Dorlf, at least in the immediate vicinity 
of the factorizing field B s , where they will be described by the state (|28|) and will therefore 
be independent of both separation and coupling range. 
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0.5 B s 1.0 1.5 2.0 0.5 B s 1.0 1.5 2.0 

B/J x B/J x 

FIG. 2. Left: The quantum discord D and the geometric discord 1% between spin pairs in the exact 
ground state of a fully connected XY spin 1/2 array of n = 50 spins with coupling anisotropy 
X = Jy/Jx = 0.5, as a function of the scaled transverse applied field B. The dotted lines depict 
the result for the mixture of aligned states (|28[) at the mean field angle cos 9 = B/J x , which is 
almost coincident with the exact result for \B\ < J x and exactly coincident at the factorizing field 
B s = y/xJx- Right: Same quantities for a cyclic chain of n = 50 spins with first neighbor XY 
couplings with the same anisotropy, for a distant pair (L = 25). The dotted lines depict again the 
results for the mixture of aligned states (|28p at the mean field angle, which now coincide with the 
exact results only at the factorizing field B s . 



As illustration, we first depict in Fig. [2] the quantum discord and the geometric discord 
of spin 1/2 pairs for XY couplings (J l J = V in the ground state of the fully connected 
array (Lipkin-type model) and of the nearest-neighbor cyclic chain, for anisotropy J y / J x = 
0.5 in the ferromagnetic type case (J x > 0). The emergence of a finite appreciable value 
of these quantities for \B\ < B c , persisting for pairs with large separation even in the case 
of first neighbor couplings, is then a direct consequence of the definite parity effect. The 
exact results for n = 50 spins were computed by direct diagonalization of H in the Lipkin 
case (where the ground state belongs to the completely symmetric representation having 
total spin S = n/2), while in the nearest neighbor chain they were obtained through the 
exact Jordan- Wigner fermionization 66], taking into account the parity effect exactly in the 



discrete Fourier transform (see for instance Refs. |43|, |61|, |62] ) . In both cases the ground state 
exhibits n/2 parity transitions as the field increases from 0, the last one at B s , although for 
the case depicted (N = 50), their effects on D or J 2 are not visible in the scale of the figure 



(they become visible for smaller n 

It is verified that at the factorizing field f|43|) . the exact results for D and l-i in both 
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0.5 B s 1.0 1.5 2.0 0.5 B s 1.0 1.5 2.0 

B/J x B/J x 



FIG. 3. The quantum discord D (left) and the geometric discord I2 (right) between spin pairs with 
separation L = 1, 2, . . . , n/2 in the exact ground state of a cyclic chain of n = 50 spins with first 
neighbor XY couplings and anisotropy x = 0-5- The results for all separations are simultaneously 
depicted. They all merge at the factorizing field B s , where they coincide with the result for the 
mixture (f28j) with cos# = B s /J x = ^fx- 

models coincide with those obtained for the mixture (128]) . i.e., with Eqs. fl3Tj) and fl32|) 
for cos9 = ^fx (Eq. (HTjl ). being then identical and the same for any pair at this point. 
Moreover, in the fully connected case the exact results for D and I2 are actually practically 
coincident with those obtained from the mixture (1281) (dotted lines) in the whole region 
\B\ < J x if 9 is the mean field angle satisfying cos 9 = B/J x (Eq. (138|) ). since the ground 
state is in this region well approximated by the definite parity states f|44|) or fl30|) . even away 
from B s . The behavior of D and I2 for B £ [0, J x ] resembles then that obtained for the 
mixture ( |28l) for 9 G [0, tt/2]. This is not the case in the chain with first neighbor couplings, 
where the agreement holds just at B s , and where D and I2 become appreciable only for 
\B\ < B c — (J x + J y )/2 (which is smaller than the mean field critical field J x but slightly 
above the factorizing field y/x)- Nonetheless, the values attained by D and I2 in the whole 
region \B\ < B s are still quite large, owing to the definite parity effect, although they exhibit 
the effects of correlations beyond the parity projected mean field description provided by 
the states f|44l) . 

We also depict in Fig. [3] results for D and I2 for all separations L = 1, . . . , n/2 of the 
spins of the pair, in the nearest neighbor coupling case (in the fully connected model they 
are obviously identical V L). It is seen that the values of D and I2 (and the same for 
^3 or other i/'s) rapidly saturate as L increases in the region \B\ < B c , reaching here a 
finite non-negligible value due to the definite parity effect, whereas for \B\ > B c they are 
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appreciable just for the first few neighbors (L = 1, 2). In the last region they can be described 



perturbatively 



20]. 



It is apparent from Figs. EH3] that the same qualitative information can be obtained 
either from the quantum discord D or the geometric discord I2, except for the type of 
maximum. That of I2 is cusp-like due to the sharp x — Y z transition in the minimizing 
measurement direction that arises as the field increases, which parallels that occurring for 
the state ( )28|) as 6 decreases (see Fig. [p. Such transition reflects the change in the type of 
pairwise correlation, and resembles that of the concurrence (which changes from antiparallel 



43|). We also mention that while the behavior of I3 (not shown) is similar 



to parallel at B s 

to I2, other // can exhibit a smoothed maximum as the transition from x to — > z in the 



measurement direction can be continuous [29]. 

We finally remark that the exact ground state pairwise concurrence in the fully connected 
case is small (of order n _1 and bounded above by 2/n), 65[ such that the entanglement 
monotones associated with D and I2 (the entanglement of formation and the squared con- 
currence) are very small in the scale of fig. EJ The same occurs with the concurrence of 

n 

largely separated pairs in the nearest neighbor case [67[, which is non-zero (but very small 
for this anisotropy and size) just in the immediate vicinity of B s 43|, |45 |. 



IX. CONCLUSIONS 

We have examined the behavior of pairwise quantum correlations in the exact ground state 
of finite ferromagnetic-type XY spin chains in a transverse field, by analyzing the quantum 
discord as well as other generalized measures of quantum correlations. We have first provided 
a brief review of the latter, which are based on general entropic forms and defined as the 
minimum information loss due to a local measurement in one of the constituents. They 
generalize the one-way information deficit and contain the geometric discord as a particular 
case, preserving at the same time the basic properties of the quantum discord like reducing 
to the (generalized) entanglement entropy in the case of pure states and vanishing just for 
classically correlated states. We have shown that all these measures indicate the presence 
of long range pairwise quantum correlations for \B\ < B c in the exact ground state of these 
chains, which arise essentially from the definite S z parity of such state and can be understood 
in terms of the model based on the mixture of aligned states (|28|) . They all reach full range 
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at the factorizing field, where they acquire a finite non-negligible constant value which is 
independent of the pair separation or coupling range and is determined solely by the coupling 
anisotropy. Such value is exactly described by the states (I2"gj) or (Hoi) , which also provide a 
quite reliable description of these correlations for all \B\ < B c for long range couplings, as 
we have seen in the case of the fully connected model. Parity effects are then seen to be of 
paramount importance for a proper description of quantum correlations in finite quantum 
systems. A final comment is that the use of simple entropic forms involving just low powers 
of the density matrix, like those underlying the geometric discord I2 and the cubic measure 
^3, enables an easier evaluation, offering at the same time an increased sensitivity of the 
optimizing measurement to changes in the type of correlation. 
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